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SCWE  APPROXIMATIONS  IH  THE  l^UMBERICAL  TREATMENT  OF  RADIATION  TRANSPORT 


THROUGH  MEDIA  WITH  FREQUENCY- DEPENDENT  ABSORPTION  COEFFICIENTS 

I.  INTRODUCTION 

investigation^  into  the  properties  of  a  numerical  treatment  of 
the  transport  of  radiation  through  slabs  of  material  is  extended.  Although 
the  integral  formulation  of  the  energy  change  equation  is  retained,  the 
method  is  made  applicable  to  a  medium  having  a  frequency- dependent  absorption 
coefficient  rather  than  the  frequency- independent  one  assumed  in  the  grey- 
atmosphere  approximation.  In  other  respects  the  formulation  is  identical 
to  that  previously  desciflfced.  The  proposed  extension  relies  on  the  possi¬ 
bility  of  making  a  valid  Approximation  to  the  absorption  coefficient  of  the 
form; 


crV,  T,  p)  =o^ 


(1) 


which  permits  a  separation  of  the  spatial  and  frequency- dependent  parts  of 
the  calculation. 

Replacing  the  exponential  integral  tables  required  in  the  previous 
formulation  are  new  tables .  These  new  functions  are  obtained  by  performing 
integrations  over  the  frequency  spectrum  of  certain  functions  of  the  absorption 
coefficient  which  single  out  the  Rosseland  or  Planck  mean  absorption  coefficient 
under  appropriate  limiting  conditions. 


II.  APPROXIMATE  FREQUENCY  DEPENDEilCE 


While  in  (1)05  (x'),  was  assumed  to  be  independent  of  0^  and  to 

depend  only  on  x',  in  fact  real  materials  display  a  strong  dependence  on 

frequency  as  typified  by  the  absorption  coefficients  of  low  -  z  materials 

2 

in  the  partially  ionized  state.  Not  only  are  the  discontinuities  at  the 
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photoelectric  edges  in  evidence  but  the  broadened  lines  may  add  a  great 
deal  of  structure  to  the  spectrum. 

As  the  temperature  and  density  of  a  material  are  changed  the  absorption 
coefficient  changes  both  in  magnitude  and  in  spectral  shape .  Over  limited 
ranges  of  temperature  and  density  It  is  sometimes  permissible  to  single  out 
the  changes  in  magnitude  of  the  spectrum  and  disregard  the  changes  in  spectral 
shape  so  that: 


(7'(V,  T,  p)  =<r^  {~>))<r2  (T,p) 

This  factorization  of  the  absorption  coefficient  will  form  the  start¬ 
ing  point  for  subsequent  treatment  of  the  frequency- dependent  case. 


III.  THE  RADIATION  FLUX 


The  expre8_Bion  for  the  radiation  flux  as  developed  in  Equation  (8) 
of  (I)  is: 

0(x)  =  2«cJ  dVj'-J^  By(x’)^  (x’)E^(j^ <7:^  dx")  dx’ 

+  1  B-j  (x')<y:^  (x')  ^  (j  <ix")  dx' 

J-00  -V' 


(2) 


where 

Bv 


o3 


Equation  (2)  takes  the  form  of  Equation  {16)  of  (l)  provided  <7^  is  not  a 
function  ofV-  Introduction  of  the  factorization  assumption,  Equation  (l), 
allows  the  frequency  integration  to  be  carried  out  inside  the  spatial 
Integral : 


dV 


(3) 
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Introducing  the  average  optical  depth, defined  In  terras  of  the  spatial 
part  of  the  absorption  coefficient 

dx*  ,  (4) 

■o 

the  flux  from  Equation  (3)  takes  the  form: 


0  (x)  =  •^  dT'  T^'  (r* )  £“2  )>  t  -  ^ 

^  1 

+  /  dz'  4  (TCtr* ),  r-  zr ) 

-'•OO 


(5) 


In  Equation  (5)  the  new  kernel  has  been  defined  as  follows: 


£0  (T,x)  = 


rOO 

I  db/Ejj  (T)^ 

— 

/  dV  BV  (T) 


(6) 


In  case  =  1,  <?2  ■  (^)  as  required  by  the  grey  approximation. 

The  essential  simplification  introduced  by  factorizing  the  absorption 
coefficient  occurs  in  the  argument  of  the  function  permitting  a  single 
optical  depth  to  enter  as  a  parameter  there .  The  function  ^  as  defined  in 
Equation  (6)  depends  on  two  parameters.  The  argument  x  is  the  optical  depth 
from  source  to  point  at  which  the  flux  is  being  evaluated.  In  addition,  T, 
the  terai^erature  at  the  source  point,  will  change  6^  if  the  frequency 
dependence  differs  strongly  from  constant  over  the  range  of  the  weighting 
function. 


IV.  HUMEg^ICAL  IICTBGRATION 

In  order  to  evaluate  the  integi’al  in  Equation  (5)  numerically,  it  may 
be  desirable  to  assume  that  the  temperature  is  constant  within  a  zone.  This 
amounts  in  Equation  (3)  to  assuming  a  spatially  constant  sovurce  through  the 
zone .  Carrying  out  the  integral  over  each  zone  and  summing  all  of  the  zone 
contributions  gives  the  results; 
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0  (Xj) 


2kC 


The  zones  have  heen  chosen  In  the  same  way  as  in  (l)  with  the  flux  being 
evalviated  at  an  interface  between  zones.  In  terms  of  the  optical  depth  "Zr 
of  Equation  (4): 


0  (Xj)  =  ^  ^i+1/2^  [^3  ^'^^^i+1/2^’  ^1  " 

^  (T(tri+2./2^>  '^.+1  ■ 


X  ^i+1/2^  ^3  ^'^^^i+1/2^'  ‘  ^+1^  ‘  ^3  ^^^^i+1/2^^ 


J;1 

+ 


where 


63  (T,x)  = 


,00 

J  dV  Bv  (T)  x) 

f  dV  Bv  (T) 


(8) 

(9) 


When <7^  is  independent  of  frequency  Equation  (8)  reduces  to  the  exact 
form  proposed  in  (l)  in  zones  removed  from  the  point  at  which  the  flux  is 
being  evaluated.  While  the  prescription  indicated  above  is  suitable  for 
optically  thin  zones  and  for  optically  thin  material,  a  more  refined 
description  of  the  temperature  between  zones  is  required  to  obtain  the  proper 
limiting  behavior  for  optically  thick  zones  and  material.  The  next  order 
of  refinement  is  to  express  the  source  function  as  a  linear  function  of 
optical  depth  in  the  neighborhood  of  the  point  at  which  the  flux  is  being 
evaluated: 


dx) 
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The  contrihution  of  the  two  zones  adjacent  to  the  Interface  at  which 
the  flux  is  evaluated  is  replaced  by  that  obtained  by  assuming  the  distribution 
of  Eqviation  (lO): 

“j  [®3  ■  ®3 

+  o-^  ^j+1  ■  ^l^^+l  ‘^j^^ 

+  er^  (ZTj  -  ^.tej  - 

Carrying  out  the  frequency  integrations  gives: 

0  (xj)  =  —  <  - Vj.i  "^j+i/s  ^3  ^"^j+i/s '  "^j+i  ■ 

■  ^3  ^^J+1,/2  '  ^  ■  ^-1^] 

^  ^+1  ■  ^-1  ^J-V2  [^3  ^'^J-l/2  '  ^J+1  ■  ■  4  ^^J-l/2  ’  ■  '^-1^] 


0(Xj)  =  2ncJ  dT^f 
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+  2 


*  ^4  (■'j  •  rj  -  25.1)  (12) 

*  <^J+l  -  ZTj)  ^3  (fj  .  ■  tj)  *  (Cj  ■  (Tj  ,  ZJ  -  1  , 


where 


(T)  H 


Jo  "^1  ^ 


,-l 


^  755 — as: 


dH-? 

dT 


+3^  (T,  0) 


1  -1 


£  (T  _  X  ^1  d/  ^ 
H  ,00  dB,, 

Jo  dr 


•'  (T,  x)  =-^ 


00  ^B-' 


dV— ?  E  (tT-  x) 

o  dT^ 


^00  dB  . 

[  -I 

Jo  dT 


Equation  (12)  Is  more  complex  than  the  corres^'r'^^ng  one  for  the  grey 

atmosphere  approximation  hut  reduces  precisely  to  it  whencr^  is  independent 

of  V.  In  that  case  6,  =  =  E_,  so  that  only  two  functions,  E_  and  Er,  are 

required.  It  Is  possible  that  the  difference  between  and  is  of 

secondary  importance.  The  Planck  mean, <7%,,  of  the  frequency  dependent  factor 

n 

is  the  leading  term  of  the  expression.  Since  the  linear  Interpolation  is 
important  only  for  optically  thick ^  zones  the  contributions  of 
small.  It  is  likely  that  either  could  be  replaced  by  ^  or  the  inter¬ 
polation  could  be  modified  to  achieve  the  replacement .  Whether  simplification 
of  these  terms  is  possible  is  to  be  determined  by  test  calculations  yet  to 
be  performed. 


In  order  to  carry  out  the  calculation  It  Is  proposed  that  tables  of 
prepared  to  be  used  instead  of  and  Ej^.  The  table 
prepai-ation,  while  a  bigger  Job  than  finding  the  Roseeland  and  Planck  means, 
can  be  computed  with  relatively  little  difficulty  and  computer  time  before 
the  time- dependent  part  of  the  transport  calculation  is  performed. 


V.  PROPERTIES  OF  THE  FUNCTIONS 

The  6- functions  entering  the  frequency- dependent  formulation  display 
properties  which  depend  on  various  means  of  the  absorption  coefficient  as  a 
function  of  frequency. 

The  emission  of  energy  from  an  optically  thin  region  will  depend  on 
the  ^  fxinctlon.  The  value  of  for  small  argument  is: 

(T,  x)  =  1/2  -  (T)  ,  (13) 


where  o-p  is  the  Planck  frequency  mean  of  the  frequency  factor of  the 
absorption  coefficient.  In  the  case  of  very  transparent  material  Equation  (8) 
for  the  flux  out  of  the  material  reduces  with  Equation  (13)  to  a  svun  of 
terms.  Each  term  represents  the  Planck  emission  of  a  zone  of  the  material 
and  contains  the  frequency  dependence  correctly  through  the  Planck  mean. 

For  regions  of  material  which  are  optically  very  thick  the  corrections 
given  in  Equation  (12)  are  applicable.  The  leading  term  and  the  only  one 
which  survives  for  very  thick  zones  is: 


f  (Xj)  . 


which  is  the  diffusion  approximation  to  the  radiation  flux  containing  the 
Rosseland  mean  of  the  frequency  factor.  The  next  term  ^  is  a  generalization 
of  the  Rosseland  mean  which  weights  even  more  strongly  the  regions  of  small 
absorption  coefficient.  However,  it  and  the  ^  terms  will  be  relatively 
small  in  the  diffusion  regime  by  virtue  of  the  exponentially  decreasing 
character  of  the  E„  and  El  functions  of  large  argument .  The  function 
contains  the  derivative  of  the  Planck  function  as  does  but  contains  E^ 
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and  has  noo^^.  appearing  outside  the  argiunent.  In  this  sense,  it  represents 
a  combination  of  optically  thin  and  thick  properties.  Its  role  is  not  yet 
understood . 


VI.  SUGGESTIONS  FOR  FUTURE  DEVELOFMBMT 

The  method  proposed  offers  the  possibility  of  treating  frequency 
dependence  by  a  numerical  treatment  with  the  equivalent  of  a  single  frequency 
group,  having  very  mild  stability  restrictions,  and  giving  reasonable  accuracy. 
It  is  attractive  as  compared  to  alternatives  which  call  for  solution  of  the 
intensity  equations  with  attendant  large  storage  requirements,  possibly  small 
time  steps  associated  with  photon  transit,  angular  groups,  and  frequency 
groups.  In  particular,  the  latter  has  appeared  to  require  a  large  prolifer¬ 
ation  of  the  calculation  in  order  to  adequately  treat  the  complexity  presented 
by  the  absorption  coefficients  of  real  materials. 

However,  the  approximation  Equation  (l)  upon  which  the  frequency  treat¬ 
ment  rests  may  require  refinement  to  treat  general  problems.  Several 
possibilities,  none  of  which  has  been  investigated  to  any  degree,  suggest 
themselves  as  being  relatively  easy  and  inexpensive  extensions  of  the  present 
method: 

1.  Divide  the  spectrum  into  coarse  groups,  in  each  of  which 
the  factorization  is  assumed  to  hold  but  with  different 
spatial  dependence. 

2.  Express  the  absorption  coefficient  as  a  sum  of  factored 
terms,  the  spatial  dependence  being  different  in  each 
term. 

3.  Expand  the  dependence  of  absorption  coefficient  on 
temperature  and  density  about  a  central  value  and  retain 
the  constant  and  linear  terms. 
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